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Abstract. We present new spectroscopic observations of 
Mkn 309, a starburst galaxy with one of the largest WR 
populations known. A highly super solar metallicity of 12 
-|- log(0/H) ^ 9.3-9.4 is derived. Using additional ob- 
jects from Guseva et al. (2000) we analyse a sample of 
five metal-rich ([0/H] > 0) WR galaxies with the main 
goal of constraining the basic properties of the massive 
star populations (IMF slope, M^p) and the star formation 
history (age, burst duration) of these objects by quantita- 
tive comparisons with evolutionary synthesis models. The 
following main results are obtained: 

— The observations are well explained by extended bursts 
of star formation with durations At ~ 4-10 Myr seen 
at ages of 7-15 Myr or a superposition of several bursts 
with age differences of ~ 4-10 Myr including a young 
(^ 5 Myr) burst. This naturally explains both the ob- 
served WR populations (including WN and WC stars) 
and the presence of red supergiants. The burst dura- 
tions, somewhat longer compared to those derived in 
other WR galaxies using the same models (Schaerer et 
al. 1999a), are plausible in view of the physical sizes 
of the observed regions and the nature and morphol- 
ogy of our objects (nuclear starbursts), and pose no 
fundamental physical problem. 

— The SEDs in the optical range are very well reproduced 
for all objects, provided the stellar light suffers from a 
smaller extinction than that of the gas (derived from 
the Balmer decrement). This confirms earlier findings 
from studies combining UV-optical data of other star- 
burst galaxies. 

— All the considered observational constraints are com- 
patible with a Salpeter IMF extending to masses M^p 
^ 40 Mq. Adopting a conservative approach we de- 
rive a lower limit of Afup^ 30 Mq for the Salpeter 
IMF. From more realistic assumptions on the metal- 
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licity and SF history we favour a lower limit M^p^ 
30-40 Mq, which is also in agreement with II/3 equiv- 
alent width measurements of metal-rich H ii regions in 
spiral galaxies indicating an upper mass cut-off of at 
least ~ 35 - 50 Mq. Steep IMF slopes (a > 3.3) are 
very unlikely. 

The uncertainties of our results are discussed. We com- 
pare our findings to other work on massive star popula- 
tions and the IMF in similar environments. We stress the 
importance of direct analysis of stellar populations com- 
pared to other indirect methods based on properties of 
ionized gas to constrain the IMF in mctal-rich starbursts. 

Key virords: Galaxies: abundances - Galaxies: evolution 
- Galaxies: ISM - Galaxies: starburst - Galaxies: stel- 
lar content - Stars: luminosity function, mass function - 
Stars: Wolf-Rayet 



1. Introduction 

Massive galaxies, the central regions of galaxies, AGNs, 
and related objects are generally expected to harbour 
metal-rich gas and stars, as indicated by numerous ob- 
servations (e.g. Hamann 1997, Henry & Worthey 1999). 

For the understanding of actively star-forming galax- 
ies (also loosely denoted as "starbursts" here) knowledge 
about massive stars — their evolution, as well as radiative, 
mechanical, and chemical feedback — and their relative 
number compared to lower mass stars is of prime impor- 
tance. However, the properties and evolution of individual 
metal-rich high mass stars, massive star populations, and 
their initial mass function (IMF) in metal-rich environ- 
ments are still poorly known. 

Although of great interest to a number of astrophysi- 
cal problems few quantitative studies of individual massive 
stars at high metallicities exist so far. See e.g. McCarthy 
et al. (1997) and Monteverde et al. (1997) for work on 
M31 and M33 stars, Najarro et al. (1997) and Figer et al. 
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(1994) for the use of WR/0 star statistics. Studies of inte- 
grated stellar cluster spectra, such as the pioneering work 
of Bica, AUoin and coworkers (Bica & AUoin 1986), have 
provided an extremely useful base for the understanding 
of stellar populations in more complex objects. However, 
templates of both young and metal-rich clusters (metal- 
licities above ^ 1.25 x solar) are missing. This clearly 
shows that studies of massive stars, both individual stars 
or integrated populations, with super-solar metallicity re- 
main largely unexplored territory. While in principle neb- 
ular lines can be used through photoionization modeling 
to constrain the stellar content of H ii regions (e.g. Garcfa- 
Vargas et al. 1995, 1997, Stasihska & Leitherer 1996) the 
inherent difficulties of optical studies at high metallicities 
(e.g. strong sensitivity of lines to unknown electron tem- 
perature) render this approach rather uncertain. Here, we 
preferentially rely on the analysis of direct stellar features. 

Of particular interest for studies of metal-rich star- 
bursts is the upper part of the IMF. More precisely, im- 
portant questions are what is the slope of the IMF and 
what are the most massive stars formed in such environ- 
ments? From work undertaken in recent years the picture 
of a "universal" IMF with a slope close to the Salpeter 
value for the mass range of interest here (M ^ 5 Mq) and 
the existence of massive stars (Mup~ 60-100 Mq) seems to 
emerge (e.g. Larson 1998, references in Gilmore & Howell 
1998, but also Scalo 1998). So far these results are found 
to hold for a diversity of objects and metallicities from so- 
lar to very metal poor environments (SMC, blue compact 
dwarf galaxies with 12 -f- log(0/H)^ 7.2). Very little is, 
however, known about the IMF at high metallicities. 

Indirect indications for a possible systematic differ- 
ence of the IMF come e.g. from studies on IR-luminous 
starburst galaxies, where generally a soft radiation field 
is observed, which, among other explanations, can be in- 
terpreted as a deficit of massive stars (cf. Goldader et al. 
1997, Luhman et al. 1998). Given the high luminosity and 
intense star-formation in these objects, the metallicity is 
likely high and could be responsible for the IMF change. 
Most recently Bresolin et al. (1999) claim a low value of 
the upper mass cut-off {M^p ~ 25 Mq) at super-solar 
metallicity from considerations of the nebular properties 
of giant extragalactic H il regions. Although not yet well 
established (cf. e.g. the recent modeling of ISO observa- 
tions by Thornley et al. 2000), these results, if true, have 
many implications. It is important, however, to note that 
the only indications for a low value of M^p in these ob- 
jects come from the properties of the gas which is ionized 
by massive stars. Many causes (dust absorption, nebular 
geometry, uncertainties in stellar atmosphere models etc.) 
could be responsible for a misinterpretation in terms of 
Mup. A direct probe of the stellar content of metal- rich 
objects is needed. A first step in this direction is done in 
the present work. 

The so-called Wolf-Rayet (WR) galaxies (cf. Conti 
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from WR stars observed in the integrated spectra indi- 
cate the presence of massive stars (Mjnitiai ^ 25 Mq), 
provide an extremely interesting opportunity to constrain 
evolutionary models for massive stars in different environ- 
ments and to probe the upper part of the IMF. These 
issues have first been explored in the important studies of 
Kunth & Joubert (1985), Arnault et al. (1989) and Vacca 
& Conti (1992). Recent progress comes from new high 
S/N observations (e.g. Schaerer et al. 1999a - hereafter 
SCK99, Guseva et al. 2000 - hereafter GITOO) and the 
use of improved evolutionary synthesis models (Cervino 
& Mas-Hesse 1994, Meynet 1995, Schaerer 1996, Schaerer 
& Vacca 1998, hereafter SV98). Reviews on the issues ad- 
dressed in studies of WR galaxies can be found for exam- 
ple in the volume of van der Hucht et al. (1999) and in 
Schaerer (1999ab). 

The sample of GITOO includes also some metal-rich 
WR galaxies. In the present work we extend this sample 
and present a detailed analysis of five objects with super 
solar metallicity. Our main goal is to determine the ba- 
sic properties of the massive star populations (IMF slope, 
Mup) and the star formation history (age, burst duration) 
in these metal-rich starbursts. 

The structure of this paper is as follows: The obser- 
vational data is described in Sect. 2. Physical conditions 
and abundances are derived in Sect. 3. The stellar content 
of our objects is analysed in Sect. 4. In Sect. 5 we derive 
constraints on the star formation properties and massive 
star populations from comparisons with evolutionary syn- 
thesis models. The interpretation in terms of the IMF is 
discussed in Sect. 6. Our main results are summarised in 
Sect. 7. 

2. Observational data 

The observational data used in the present paper consist 
of four high metallicity {^-^ solar and above) objects taken 
from the large sample of Wolf-Rayet galaxies of GITOO 
(Mkn 589, Mkn 710=NGC 3049, Mkn 1199, Mkn 1259) 
, supplemented by the data on NGC 3049 by Schaerer et 
al. (1999a), and new observations of Mkn 309. The general 
properties of the objects are given in Table |l|. Distances 
and derived metallicities are included in Table |^ below. 
Except for Mkn 710, a bright H ii region in a barred spiral 
galaxy, the star forming regions discusses here are fairly 
compact central regions of galaxies ("nuclear starbursts"). 
Note that this is in contrast with the remaining WR galax- 
ies of GITOO, comprising mostly metal-poor blue compact 
dwarf galaxies. 

Spectroscopic observations of Mkn 309 were made with 
the Kitt Peak 4m Mayall telescope during a clear night on 
18 June 1999 with the Ritchey-Chretien RC2 spectrograph 
used in conjunction with the T2KB 2048x2048 CCD de- 
tector. We use a 2"x300" sht with the KPC-lOA grating 
(316 lines mm~^) in first order, with a GG 375 order sep- 
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Galaxy q(1950) (?(1950) tub Type 

Mkn 589 02 11 08.7 +03 52 08 14.51 S? 

Mkn 1199 07 20 28.5 +33 32 24 13.7 Sc 

Mkn 710 09 52 10.2 +09 30 32 13.04 SB(rs)ab Sbrst 

Mkn 1259 10 36 03.0 -06 54 37 13.5 SO pec? 

Mkn 309 22 50 10.0 +24 27 52 15.44 Sa Sy2 




Wavelength (A) 

Fig. 1. Observed spectrum of Mkn 309. Some heavy element absorption features, He I and WR lines are marked. All 
He I lines in blue part of the spectrum are in absorption while He i A6678 and A7065 lines are in emission. 



ination for wavelengths blueward of 7200A which is the 
wavelength region of interest here. The above instrumen- 
tal set-up gave a spatial scale along the slit of 0.69 arcsec 
pixel" "'^j a scale perpendicular to the slit of 2.7A pixel" ^, 
a spectral range of 3500-7300A and a spectral resolution 
of BA. The seeing was 1". The total exposure time was 
60 minutes and was broken up into 2 subexposures. The 
slit was oriented close to the parallactic angle, so that no 
correction for atmospheric differential refraction was nec- 
essary. 



The Kitt Peak IRS spectroscopic standard star Feige 
34 was observed for flux calibration with a 6" wide slit. 
Spectra of He-Ne-Ar comparison arcs were obtained be- 
fore and after each observation to calibrate the wavelength 
scale. 

The two-dimensional spectra were bias subtracted and 
flat-field corrected. We then use the IRAfQ software to 
perform wavelength calibration, correction for distortion 

^ IRAF is distributed by National Optical Astronomical Ob- 
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and tilt for each frame and extraction of one-dimensional 
spectra. The extracted spectra from each frame were coad- 
ded and cosmic rays were removed manually. The sensitiv- 
ity curve was obtained from the observed spectral energy 
distribution (SED) of Feige 34 with a high-order polyno- 
mial. The intensity distribution of the continuum and of 
the emission lines along the slit shows a strong maximum 
with FWHM = 2". We extract a one-dimensional spec- 
trum of the brightest part of the galaxy in a 2"x4'.'2 aper- 
ture. The signal-to-noise ratio for this spectrum is ^ 100. 
The spectrum for Mkn 309, uncorrected for extinction, is 
shown in Figure 

We detect several Wolf-Rayet emission lines. The He 
I lines in the blue part of the spectrum are in absorption 
suggesting the presence of early B stars in the star- forming 
region. We also mark several absorption features of heavy 
elements, such as Ca ii A3933, Na i A5889, 5895, the TiO 
7' band at A6250 and the Fe i Ba 11 + Ca i A6495 blend. 
Although weak, the latter feature is not an artifact. It is 
seen in the spectra of all studied galaxies (Fig. 2) because 
of their very high signal-to-noise ratio. Although Bica & 
Alloin (1986) note that the latter blend (their window 
#59) includes also CN and TiO features, the narrowness 
of the blend suggests that contribution of molecular bands 
is rather small. 

The observed line intensities have been corrected for 
interstellar extinction using the reddening law by Whit- 
ford (1958). We decided not to correct line intensities sep- 
arately for the Galactic and intrinsic extinctions, partly, 
because of the large discrepancy between the determina- 
tions of Ab = 0.649 by Schlegel et al. (1998) and Ab 
= 0.190 by Burstein & Heiles (1982). Instead we derive 
extinction directly from the observed hydrogen line in- 
tensities which have been also corrected for underlying 
stellar absorption, with the equivalent width for absorp- 
tion hydrogen lines derived self-consistently together with 
the extinction coefficient from the observed intensities of 
aU hydrogen lines. The redshift z = 0.04174±0.00009 de- 
rived from the observed wavelengths of 16 emission lines 
is close to z = 0.04215 in the NED. For the latter helio- 
centric redshift, 1" = 817 pc, assuming Hq = 75 km s^^ 
Mpc"-^. We show the observed F{\)/F(H(3) and extinc- 
tion and absorption-corrected /(A) / /(H/3) line intensities 
for Mkn 309 in Table ||, along with the extinction coeffi- 
cient and the equivalent width of the hydrogen absorption 
lines, the observed flux and equivalent width of the H/3 
emission line. The errors of the line intensities include the 
errors in the placement of the continuum and the errors 
introduced by Gaussian fltting of the line profiles. The 
indicated uncertainties represent formal errors. Note the 
small intensities of the [O i] and [S 11] emission lines sug- 
gesting that Mkn 309 is a starburst galaxy, contrary to its 
classification in the NED as Sy2 galaxy. 
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Table 2. Emission line intensities in Mkn 309 



Ion 


F( X) / FCHB) 




Q797 fo ttI 


n 40+0 09 


(SI +0 0"^ 


4101 HA 


n 1 0+0 09 


9^+0 0^ 




n "^1 +0 09 


4^+0 0"^ 


4:tJUtJ IVV lilj 


03+0 01 


03+0 01 


4861 H/3 


1.00±0.02 


1.00±0.03 


4959 [O III] 


0.15±0.02 


0.13±0.02 


5007 III 


0.38±0.02 


0.34+0.02 


5056 Si TT 


0.05±0.02 


0.04+0.02 


5755 fN III 


0.05±0.01 


0.03±0.01 


5876 He i 


0.02+0.01 


0.02+0.01 


6300 [O i] 


0.05±0.01 


0.03±0.01 


6563 Hq 


5.41+0.10 


3.04+0.06 


6583 [N 11] 


3.83±0.07 


2.12+0.04 


6678 He i 


0.02+0.01 


O.OliO.Ol 


6717 [S 11] 


0.21+0.01 


O.lliO.Ol 


6731 [S 11] 


0.27±0.01 


0.15±0.01 


7065 He i 


0.02+0.01 


O.OliO.Ol 


7135 [Ar III] 


0.05±0.01 


0.02+0.01 


C(H/3) dex 


0.68±0.02 


i^(H/3)" 


1.74+0.11 


EW(}if3) A 


11+1 


EW{abs) A 


O.SiO.l 


"in units of 10" 


-14 —1 

ergs s cm 


-2 



3. Physical conditions and heavy element 
abundances in the H II region 

From the observed spectrum we shall now derive the phys- 
ical conditions and metal abundances in Mkn 309. For 
the remaining objects the reader is referred to GITOO. 
Given the weakness or absence of [O ill] A4363 in Mkn 
309 and other objects included later for comparison, the 
abundance determinations are affected by relatively large 
uncertainties. These are briefly discussed below. 

3.1. Mkn 309 

To determine heavy element abundances, we adopt, fol- 
lowing Izotov et al. (1994, 1997), a two-zone photoionized 
H II region model: a high-ionization zone with tempera- 
ture Te(0 III), where O ill and Ne III lines originate, and 
a low-ionization zone with temperature Te(0 li), where O 
II, N II and S 11 lines originate. As for the Ar ill lines, they 
originate in the intermediate zone between the high and 
low-ionization regions (Garnett 1992). 

A weak auroral [O ill] A4363 emission line is detected 
with the peak intensity above the 3(T level. However, the 
intensity of this line is rather uncertain due to underly- 
ing H7 stellar absorption and uncertainties in the place- 
ment of the continuum because of numerous absorption 
features. The temperature Te(0 ill) calculated using the 
[O III] A4363/(A4959-KA5007) ratio is unreasonably high 
~ 47000 K. The electron temperature derived from the [N 
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is still significantly larger than the electron temperature 
Te(0 III) derived from the empirical relation by Edmunds 
& Pagel (1984) and can also be subject to the enhance- 
ment mechanisms discussed above. A similar re(N ii) has 
been derived by Osterbrock & Cohen (1982). Adopting 
Te(N ii) = 9900K results in an oxygen abundance 12 -f 
log 0/H = 7.85. However, the N/0 abundance ratio in 
this case is 30 - 40 times larger than the value expected 
for such a low oxygen abundance, which seems unreason- 
able. Therefore, we do not use the [O ill] A4363 and [N 
ll] A5755 auroral lines for the determination of electron 
temperature. Instead, the empirical method by Edmunds 
& Pagel (1984) is used to derive Te{0 III) from the inten- 
sities of nebular [O ll] A3727 and [O ill] A4959 + A5007 
lines. The electron temperature for [O ii] is derived using 
the results of the H II photoionization models of Stasihska 
(1990) and the electron temperature for [Ar ill] is derived 
using the prescriptions by Garnett (1992). The electron 
number density is calculated from the [S ii] A6717/A6731 
line intensity ratio. Then Tg and are used to derive 
element abundances. 

In Table || we show the derived electron tempera- 
tures Te for different ions and the electron number density 
Ne- The electron temperatures are particularly low, much 
lower than those derived from the [O ill] and [N ll] line 
intensity ratios. We find Ne = 1110 cm~'^. The resulting 
abundances for heavy elements are shown in Table ^ The 
oxygen abundance 12 -I- log 0/H = 9.41 is ~ 3 times 
larger than the solar value. In Table js] we also show the 
oxygen abundance derived from empirical methods using 
[O ii] + [O III] intensities (Edmunds & Pagel 1984) and 
the [N ii] A6583/Ha line intensity ratio (van Zee et al. 
1998). All methods give nearly the same value of oxygen 
abundance. The nitrogen-to-oxygen abundance ratio is 
4 times larger than the solar value and within the errors 
is consistent with the secondary nitrogen production ex- 
pected at the high metallicity of Mkn 309. 

3.2. Uncertainties 

Although various empirical methods concur to indicate a 
similar O /H abundance for Mkn 309 (cf. above) , metallic- 
ity determinations in the absence of direct electron tem- 
perature measurements remain uncertain. This is the case 
for most of the objects in this study. The quoted errors 
indicate only formal uncertainties. 

McGaugh (1994) estimates an uncertainty of ^ 0.1 dex 
on 0/H at high abundances. As shown by Stasihska (1998) 
calibrations of the strong line method are also affected by 
age and geometrical effects. Errors are difficult to estimate 
in this case; from theoretical modeling an uncertainty of 
0.1-0.2 dex may not be unreasonable (Stasihska 1998). 

In GITOO and the present paper we adopt the van Zee 
et al. (1999) calibration with typical errors of 0.2 dex and 
which is formally only valid to 12 + log(0/H)< 9.1. This 



Table 3. Heavy element abundances in Mkn 309 



Parameter 



Value 



n{0 iii)(K) 
mo II) (K) 
re(Ar iii)(K) 
iV.(S n)(cm-3) 

O+/H+(xl0^) 
O++/H+(xl0=') 
O/H(xl0=*) 
12 + log(0/H) 

12 + log(0/H) (Edmunds & Pagel) 
12 + log(0/H) (van Zee et al.) 

N+/H+(xlO^) 

ICF(N) 

log(N/0) 

Ar++/H+(xlO^) 

ICF(Ar) 

log(Ar/0) 



3830±520 
6110±520 
4970±520 
1110±350 

0.38±0.21 
2.20±2.29 
2.58±2.30 
9.41±0.39 
9.32±0.20 
9.32±0.20 

0.20±0.10 
6.85 

-0.29±0.64 

1.68±0.75 
2.25 

-2.84±0.41 




-| , , 1 r- 

Mrk 309 

Mrk 589 



TiO 7' Fel + Ball + Cal 



Mrk 1259 



TiO 7 



6000 



6500 7000 
wavelength [A] 



7500 



Fig. 2. Spectral region showing TiO bands. From top to 
bottom: Mkn 309, 589, 710, 1199, 1259. The narrow Fe i 
-t- Ba II -I- Ca I A6495 blend is also marked. 



mentioned above and affects other calibrations in the lit- 
erature (e.g. Edmunds & Pagel 1984). The van Zee et al. 
calibration is based on the intensities of [N ll] lines. Our 
metallicities could be overestimated if nitrogen is enriched 
in nuclear regions of galaxies (e.g. Coziol et al. 1999). 
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Table 4. Observed parameters of WR lines 



Parameter 


Value" 




2.20±0.17 


F(N III A4512) 


0.25±0.03 


FiSi III A4565) 


0.06±0.02 


F(N V A4612) 


0.14±0.02 


F(N III A4640) 


0.25±0.02 


F{C IV A4658) 


0.22±0.02 


F(He II A4686) 


0.23±0.02 


F(C III A5696) 


0.10±0.02 


F{C IV A5808) 


0.15±0.02 


£^■(4650)" 


4.09±0.23 A 


EW{C IV A5808) 


1.19±0.20 A 


"in units of 10 ergs s ^cm 
'measured in a 6" x 6" aperture 



''includes the three WR hues from 4640 - 4686 A. 



4. Stellar features and analysis of the stellar 
populations 

The quality of the spectra allow the detection of several 
stellar emission and absorption features the most promi- 
nent ones being the WR emission lines, H and He absorp- 
tion lines, TiO, Ca ii, Na i and Fe i absorption features. 

The broad emission lines attributed to WR stars are 
discussed below and will later be used extensively for 
quantitative modeling. The high order Balmer absorption 
lines as well as He i absorption lines, often seen in H ii 
regions, are mostly due to OBA stars. Recent synthesis 
models including these lines are given by Gonzalez Del- 
gado et al. (1999; cf. also Olofsson 1995). 

4-.1. Features from late type stars 

Several broad absorption features, mostly due to TiO 
bands, are observed around ~ 6250 and ~ 7200 A in all 
our objects (see Fig. ||) . These absorption features as well 
as some narrow absorption lines of atoms and ions, com- 
monly observed in integrated spectra of clusters including 
young stars (cf. Bica & AUoin 1986), indicate the pres- 
ence of red giants or supergiants in our regions. However, 
some contribution of interstellar Ca ii and Na i absorption 
might be present. Quantitative measurements of TiO ab- 
sorption features for a large sample of blue compact dwarf 
galaxies have been provided by Guseva et al. (1998)[[ The 
two absorption features found in our objects are observed 
in giant and supergiant stars of spectral types ~ K3 and 
later (e.g. Silva & Cornell 1992). Despite the presence of 
nebular lines there is no indication of the broad TiO 7' 
band at ~ 6540 - 7050 A seen in M types. The observed 

^ Using the TiO features, stellar metallicities have been ob- 
tained for 36 star-forming regions. At least for the lowest metal- 
licities, their results indicate that the metallicities of stars are 
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TiO features are thus compatible with giants or super- 
giants of spectral type K3 or later K types. Presently 
we are not able to firmly distinguish between a young (Jj 
7-10 Myr) supergiant population or an older giant popula- 
tion (possibly a "nuclear" or bulge like population). How- 
ever, the equivalent widths of heavy element absorption 
features are rather consistent with young stellar popula- 
tion. While strength of TiO bands is insensitive to the age 
of stellar population (Bica & AUoin 1986), the equivalent 
widths of Ca 11 A3933 (1.5 A) and Na i A5889, 5895 (1 A) 
lines are smaller than those for old stellar population and 
are consistent with values expected for the young stellar 
population (Bica & Alloin 1986). Therefore, we strongly 
favour the first interpretation (population of red super- 
giants), which, as will be shown later, allows a consistent 
explanation of all major observational constraints in the 
framework of a simple star formation history dominated 
by recent activity. Unless otherwise stated this working 
hypothesis is thus adopted in the remainder of this work. 
A more firm distinction should be possible using obser- 
vations covering a more extended wavelength range and 
other stellar signatures. 

Limited quantitative information can be obtained from 
the observed TiO bands. First contamination by nebular 
emission lines and sky lines render exact measurements of 
some features difficult. Furthermore, quantitative models 
of the TiO features in young populations are not available 
yet. Although promising for future work (cf. Schiavon et 
al. 2000), no detailed quantitative analysis of the TiO fea- 
tures will thus be undertaken here. The only constraint 
retained below (Sect. 5), is that we conclude that the RSG 
features indicate the presence of stars with ages 7-10 
Myr in the observed region (e.g. Origlia et al. 1999). 

In passing we note the similarity of some of our find- 
ings with earlier studies on the link between starbursts and 
AGNs, which have revealed the simultaneous presence of 
the Ca II triplet originating from red giants or supergiant, 
H ii-region like emission lines, and partly also WR stars 
(e.g. Terlevich et al. 1990, 1996, Pastorizaet al. 1993). The 
main aim of our study being different from these studies, 
we leave a discussion of similarities and differences to fur- 
ther investigations. 

4-. 2. Broad emission lines from WR stars - analysing the 
WR population 

The objects included in our study have been selected for 
the presence of broad emission lines attributed to WR 
stars. The following lines have been identified in our ob- 
jects: N III A4512, Si III A45650, the so-called blue bump: 
N V A4612, N III A4640, C iv A4658, He 11 A4686; C iii 

^ On further inspection of the data of GITOO and additional 
observations, it appears that the origin of this line (stellar or 
nebular) may not be fully established. This will be addressed 

:„ _ -C—x Ul; x: 
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A5696, and the red bump C iv A5808 (see SCK99, GITOO). 
The WR line intensities and equivalent widths in Mkn 589, 
Mkn 710, Mkn 1199 and Mkn 1259 are taken from GITOO. 
Additional data on Mkn 710 from SCK99 are also included 
our analysis. The new data on Mkn 309 is presented be- 
low. The intensities of the blue WR bump in those and 
other objects to be discussed later are also given in Table 
. For objects in common with earlier studies, our intensi- 
ties (/(WR)//(H/3)) and EW(H(3) are in good agreement 
with the data of Mas-Hesse & Kunth (1996) and SCK99 
for Mkn 710, and with that of Ohyama et al. (1997) for 
Mkn 1259. Differences with the data of Mkn 710 from 
Vacca & Conti (1992) were discussed in SCK99; their ori- 
gin remains unclear. 

We shall now discuss the observed WR lines and esti- 
mate the number of WR stars of various subtypes in Mkn 
309. Such an analysis is given by SCK99 and GITOO for 
the other objects. 

Several broad emission lines (FWHM - 15 - 50 A) 
have been detected in the spectrum of Mkn 309 suggesting 
the presence of Wolf-Rayet (WR) stars of several types. 
The observed fluxes of these lines measured in the 2"x4'.'2 
aperture are shown in Table ^. The H/3 flux given there was 
measured in a 6"x6" aperture using another spectrum. 
Comparison with the narrow slit spectrum shows that the 
large aperture contains most of the H/? emission which is 
more spatially extended than the continuum emission. 

These fluxes have then been corrected for extinction, 
and line luminosities were calculated adopting the dis- 
tance D = 168.6 Mpc derived from z=0.04215 (NED) 
assuming Ho^75 km s~^ Mpc^^. The luminosities have 
been calculated for two cases: I) the extinction coefficient 
C(H/3) is equal to 0.68 as derived from the Balmer decre- 
ment for both gas and stars and II) C(H/?) is equal to 
0.68 for the gas and 0.15 for stars (Sect. 5). The results 
are shown in Table |5[ The H/3 flux and equivalent width 
agree well with the data of Osterbrock & Cohen (1982, 
hereafter OC82). The intensity of their WR bump is ^ 
30 % larger than our measurement including all WR lines 
from 4512 to 4686 A. Such differences could have various 
origins (different observational setups etc.) and are not 
uncommon (see e.g. comparison in SCK99). They do not 
affect our conclusions on Mkn 309. 

The strongest lines are those of the most frequently 
observed blue "WR bump" ('^ 4650 A) which we identify 
as a composition of the N v A4603,4619 doublet denoted 
here as N v A4612, the Niii A4634-40 blend, the Ciii/iv 
A4650-58 blend, and He ii A4686. Broad N ill A4512 is also 
detected as previously noted in several objects by GITOO. 
We confirm the presence of N v A4612, C ill A5696 and 
C IV A5808 which were suspected by OC82. 

The observed WR features clearly indicate a mixed 
population of WN and WC stars. What can be said about 
their subtypes ? The low N v A4612/N ill A4640 ratio 
indicates a late WN subtype; the observed value corre- 



with subtypes WN7-8 (Smith et al. 1996). C iii A5696 
emission is only strong in late WC stars. The observed 
C III A5696/C IV A5808 ratio, the primary classification 
criterion, corresponds to WC7-8 stars (SV98, Crowther et 
al. 1998). Taking into account the contribution from the 
detected WNL stars to C IV A5808 changes the subtype 
to WC8 according to the Crowther et al. classification. 
The width of C ill A5696 (FWHM ~ 19 A) is, however, 
smaller than that of C iv A5808 (49 A). This could in- 
dicate a mixture of late and some early WC stars. Our 
improved observations confirm the presence of WN and 
WC stars which were suspected by OC82 and re-enforce 
their proposed subtype distribution. 

We now estimate quantitatively the number of WR 
stars using the calibration of WR line luminosities of SV98 
from Galactic and LMC stars. The resulting numbers of 
WR stars and their relative numbers are shown in Table 
^ for the two choices of extinction coefficient mentioned 
above. Since He ii A4686 is predominantly emitted by WN 
stars and no nebular emission is present given the low 
excitation we use this line to count the number of WNL 
stars. Depending on the adopted extinction ~ 7000-26000 
stars are found. The number of WC stars can be derived 
from C IV A5808 and/or C ill A5696. For the former line 
the contribution of emission from WNL stars 10-15 %) 
has to be subtracted. Adopting the mean spectral type 
to be WC7, we find - 4000 - 16000 WCL. Strong vari- 
ations of the line luminosities of individual WC stars on 
subtype are found (SV98); adopting WC8 as mean spec- 
tral type would considerably increase the deduced num- 
ber of WCL stars. Assuming that both WCL and WCE 
stars are present (as indicated by the line widths discussed 
above), we can alternatively derive their numbers from 
C IV A5808 and C ill A5696 respectively. The total number 
of WCL+WCE stars derived in this way is in reasonable 
agreement with the previous estimate. In all cases the re- 
sulting N{WC) /iV(WN) ratio is found to be between ~ 0.4 
to 1.1 (Table ||). A metallicity dependence of the average 
line luminosities used here can a priori not be excluded. 
However, Galactic and LMC WR stars do not show sig- 
nificant differences, except for N ill A4640 (cf. SV98). In 
addition, since the Galactic objects used by SV98 includes 
a good fraction of WR stars located inside the solar ra- 
dius at metallicities typically ~ 0.15 dex above solar (cf. 
Maeder & Meynet 1994) our estimated WR populations 
and the detailed comparisons shown below should not be 
affected by possible variations of the properties of individ- 
ual WR stars with metallicity. 

What can be said about the remaining WR lines ? 
Using again the line luminosities from SV98, different es- 
timates of WC stars yield predictions for the C IV A4658 
line in good agreement with the observed one (differences 
from -20 to -|-50 %). The observed N III A4640 emission 
from WNL stars is found to exceed the expected value by a 
factor of ~ 2 based on the average Galactic value of N iii 
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Table 5. Derived parameters for WR lines and stellar populations of Mkn 309 



Parameter 


r 


iF 


D 

L^iN 111 A4512) 
L^lsi III A4565) 
L=(N V A4612) 
L'=(N III A4640) 
L^iC IV A4658) 
L'=(Hc II A4686) 
^^=(0 III A5696) 
L=(C IV A5808) 


168.6 Mpc 
(3.59±0.28)xl0''^ 
(4.75±0.49)xl0*' 
(1.05±0.42)xl0*" 
(2.44±0.43)xl0*° 
(4.43±0.39)xl0*" 
(3.82±0.36)xl0''" 
(4.06±0.35) x 10"*" 
(1.29±0.21)xl0*" 
(1.74±0.29)xl0*° 


168.6 Mpc 
(3.59±0.28)xl0''^ 
(1.25±0.13)xl0*'' 
(2.82±1.14)xl0^'' 
(6.65±1.18)xl0^^ 
(1.21±0.11)xl0''° 
(1.05±0.10)xl0'*° 
(1.13±0.10)xl0*° 
(4.71±0.77)xl0^^ 
(6.52±1.08)xl0^^ 


iV y\J 1 V ) 


Doyuuztoiuu 


Doyuuztoiuu 


N{WNL) 

iV(WC7) from uncorrected C IV A5808 
Af(WC7) from corrected C IV A5808 for WNL 
Ar(WC7) from C iii A5696 
iV(WC7)/iV(WNL) 


25800±2300 
12400±2100 
10600±1800 
16000±2600 
0.42 - 0.63 


7060±630 
4650±770 
4150±690 
5810±950 
0.59 - 0.82 


N(WC8) from C ill A5696 

Af(WC4) from C iv A5808 corrected for WNL+WCL 
Ar(WC4 + WC8)/Ar(WNL) 


18800±3100 
2780±460 
0.85±0.14 


6840±1120 
1140±190 
1.13±0.19 


N (WNL) /N {07V) 
nIwR)/N{07Y) 


0.39±0.05 
0.59 - 0.71 


O.lliO.Ol 
0.19 - 0.22 



"C(H/3) — 0.68 is assumed for correction of all lines. 

''C(H/3) = 0.68 and 0.15 are assumed for gaseous and stellar emission, respectively. 
°in ergs s~^. 



Table 6. Properties of related WR galaxies discussed in text (ordered by decreasing L(WR)). 



Name 


Metallicity 

[12 + log(0/H)] 


Distance 

[Mpc] 


/(WR)//(H/3) 
intensity 


L(WR) 

[org s"^] 


N{WR) 


References 


Htgh L(WR) objects 






2.9 xlO^^ 






IRAS 01003-2238 




470 


0.61 


9. xlO" 


Armus et al. (1988) 


Mkn 309 


9.3 


168.6 


0.34 


1.2 xlO*^ 


(1.2-4.) xlO* 


present work 


NGC 1614 


8.6 


64 


0.08 


7.6 xlO*° 


2.4 xlO* 


Armus et al. (1989) 














Vacca & Conti (1992) 


Mkn 1259 


8.95 


28.7 


0.27 


5.7 xlO**" 


2.0 xlO^ 


GITOO 


Mkn 1199 


9.13 


54.1 


0.38 


5.9 xlO^° 


1.9 xlO" 


GITOO 


Mkn 477 




152 


~ 0.1 


5.4 xlO*° 


1.7 xlO* 


Heckman et al. (1997) 


F08208+2816 


8.15 


189 


0.08 


3.7 xlO^° 


2.4 xlO"* 


Huang et al. (1998) 


Mkn 702 


8.60 


212.1 


0.13 


3.2 xlO^° 


1.2 xlO* 


GITOO 


Mkn 589 


8.99 


45.9 


0.24 


2.1 xlO*° 


8. xlO^ 


GITOO 


Objects with large 


I(WR)/I(HI3) 






4.6 xlO^" 






Mkn 710 


9.03 


20 


~ 0.3 


~ 800 


SCK99, GITOO 


Mkn 1063 


8.46 


20 


0.20 


1.6 xlO^^ 


~ 600 


GITOO 


Mkn 178 


7.82 


2.96 


0.39 


1.6 xlO^'' 


5 


GITOO 


NGC 5128/13 


~ 9. 


5.5 


1.0-1.8 






D'Odorico & Rosa (1982) 



a metallicity dependence (Smith et al. 1996). The excess 
could thus be due to larger N iii A4640/He ii A4686 emis- 
sion in WNL stars at high metallicities (cf. also Schmutz 
& Vacca 1999). Using the derived WR number popula- 



tion and the data of SV98, weak emission - below our 
detection limit - is expected for other WR lines. This is 
in agreement with the observed spectrum. Similar quan- 
titative comparisons cannot yet be done for N III A4512 
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and Si ill A4565 (see GITOO for detections of these lines 
in other objects). We conclude that all broad emission 
lines detections and non-detections are compatible with 
the WNL and WC populations indicated in Table |[ 




8.5 

ia + log(0/H) 



9.5 



Fig. 3. Observed /(WR)//(H/3) ratio for WR galaxies 
from GITOO (filled dots) and WR galaxies from Table 
^ (open squares and circles) as a function of mctallic- 
ity. The low L(WR) object Mkn 178 is marked separately 
(open triangle). Typical uncertainties are shown by the 
errorbars. The formal uncertainties are < 5 percent for 
the line intensities from GITOO. For objects with 12 -|- 
log(0/H) < 8.3 (measured from [O ill] A4363) the uncer- 
tainties on the metallicity are < 0.1 dex. Larger uncer- 
tainties are expected for the remaining objects (cf. Sect. 
I). NGC 5128/13 with a very uncertain /(WR)//(H/3) 
measurement is not shown. Circles indicate objects of un- 
known metallicity, placed arbitrarily at solar metallicity. 
The /(WR)//(H/3) ratios are uncorrected ior aperture ef- 
fects (cf. GITOO) and possible differential extinction be- 
tween gas and stars. The large /( WR) / /(H/3) ratio of Mkn 
178 is likely due to an exceptionally large aperture effect 
for which GITOO estimate a reduction of a factor 5 or due 
to small number statistics. 



4-3. WR and O star populations in Mkn 309 and related 
objects 

The ionizing flux of Mkn 309 corresponds to a number 
A^(O7V)w65900 of equivalent 07V stars assuming Case B 
recombination, Tg = 5000 K, and rig = 100 cm^'^. Thus, 
the number ratio of WR stars to equivalent 07V stars is 

Ar/TirT>\ / '\Tfr\'7\T\ n cn n fri ; „„„ T ^'^^„^^l„l^^ ] n 1 n 



- 0.22 in case II. Here iV(WR) = A^(WNL-HWC4-hWC8) 
or iV(WR) = iV(WNL-hWC4+WC7). The high Af(WR) / 
A^(07V) value (assuming identical extinction for gas and 
stars) is more than a factor of 2 larger than the observed 
iV(WR)/iV(0) value in Galactic regions inside the solar 
radius and in M31, which are representative for environ- 
ments with constant star formation at solar and super- 
solar metallicity (see Maeder & Meynet 1994 and refer- 
ences therein). As for other WR galaxies (cf. e.g. Vacca 
& Conti 1992, GITOO) such 7V(WR)/7V(0) ratios gener- 
ally indicate intense star forming events taking place over 
relatively short timescales (e.g. review of Schaerer 1999a). 
Estimates of the SF timescales for the objects of our study 
and allowing also for the observed lower extinction for the 
stars will be derived below (Sect. 5). 

Among aU WR galaxies with 12 + log(0/H) > 8.3, the 
largest iV(WC)/A^(WN) ratio (- 0.4 - 0.8) is derived in 
Mkn 309 (cf. SCK99, GITOO). Although lower, the upper 
value comes close to iV(WC)/7V(WN) ~ 0.9 - 1 observed 
at high metallicities in the Local Group (Massey & John- 
son 1998). As shown below, the WN and WC populations 
in all metal-rich WR galaxies studied here can be under- 
stood in terms of current stellar evolution models. 

To compare the WR populations of Mkn 309 with 
those in other related objects we have compiled in Ta- 
ble ^ the properties of WR galaxies with large WR bump 
luminosities and large /(WR)/H/3 line intensities. The 
compilation is based on the catalogue of Schaerer et al. 
(1999b) and the data discussed in Schaerer (1999a). In 
absolute terms, the WR population in Mkn 309 is among 
the largest known; presently it is only surpassed by the 
ultra-luminous infrared galaxy IRAS 01003-2238 (=IRAS 
F01004-2237, Armus et al. 1988, Veilleux et al. 1999). 
Other objects with more than - 8000-10000 WR starf^ 
are NGC 1614, Mkn 477, Mkn 1259, Mkn 1199, Mkn 702 
and Mkn 589, and F08208-f2816 whose main properties 
are summarised in Table ^ 

Interestingly, at least the two most "WR-luminous" 
objects (IRAS 01003-2238 and Mkn 309) also stand out 
by very large /(WR) //(H/3) ratios^ as shown in Fig. |[ 
The intensities (see Table ^ clearly exceed the typical 
values of <, 0.1-0.2 (cf. also Schaerer 1999a). Other ob- 
jects with strong WR bump intensities but much smaller 
absolute populations of WR stars are the H ii region # 
13 in NGC 5128 (eeCcu A), Mkn 178, Mkn 710, and Mkn 
1063 also hsted in Table | Although few such WR regions 
are currently known at the high metallicity end, it seems 



^ WR numbers based on determinations in the original pa- 
pers or objects with L(WRbump) 2.x 10*" erg s~^ using the 
cahbration of Smith (1991) for individual WN7 stars. The hst 
is likely incomplete due to missing data on WR line fluxes. 
Marginal cases not included here are NGC 6764, Mkn 315, He 
2-10 (references in Schaerer et al. 1999b). See GITOO (Fig. 3) 
for additional WR-bump luminosities. 

7(WR) stands for the intensity of blue bump including 
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that the largest WR intensities are found in massive and 
distant objects. If true, this is quite surprising. Indeed for 
similar mctallicity and IMF one would statistically expect 
reduced WR strengths in more distant objects where the 
spectroscopic observations encompass a much larger area 
containing more likely a mixture of several populations 
with different ages. Observational biases may, however, 
also be present (see Kunth & Joubert 1985). New obser- 
vations of nearby (low L) metal-rich WR regions will be 
needed to confirm this trend. In any case, it appears that 
despite the likely growing importance of the "age spread 
effect" the increase of the WR population with mctallic- 
ity (cf. Maeder & Meynet 1994) is sufficient to lead to the 
continuing increase of the /(WR)//(II/3) envelope with 
metallicity first shown by Arnault et al. (1989). 

5. Model comparison 

To constrain the basic properties of the massive star pop- 
ulations (IMF slope, Afup) and the star formation history 
(age, burst duration) of our objects we now undertake a 
detailed comparison between the observations and evolu- 
tionary synthesis models. 

5.1. Procedure 

The following main observational constraints are used: 

1. H/3 and Ha equivalent widths. The former is used as 
a primary age indicator; once EW{Il(3) is reproduced 
^^^^^(Ila) serves as an independent test for the pre- 
dicted spectral energy distribution (SED) in the red. 

2. Nebular line intensities. F (Ra) / F (HP) determines the 
extinction of the gas. The use of other line intensities 
requires detailed photoionization modeling which is be- 
yond the scope of this paper. 

3. Intensities and equivalent widths of the main WR fea- 
tures. The blue bump (henceforth referred to as 4650 
bump) and C iv A5808 (red bump) serve as main con- 
straints on the WR population. To avoid uncertain- 
ties in deblending individual contributions of the blue 
bump we prefer to use measurements for the entire 
bump. In contrast to the spectra of metallicity objects 
our low excitation spectra show no evident contami- 
nation from nebular lines (e.g. [Fe ill] A 4658, nebular 
He II A4686). 

To potentially disentangle between various effects (un- 
derlying "non-ionizing" population, loss of photons, 
differential extinction between gas and stars) it is im- 
portant to use both equivalent widths and relative 
WR/H/3 intensities (cf. SCK99). 

4. TiO bands at - 6250 and 7200 A. Assuming that they 
originate from a red supergiant population (rather 
than giants; cf. Sect. 4) this indicates the presence of 
cool stars from a population with ages ^ 7-10 Myr. 

5. SEDs available over the full range of ~ 3500-7500 A 



WR populations, which are responsible for the contin- 
uum flux. 

For the model comparisons we use calculations based 
on the evolutionary synthesis code of SV98, which in par- 
ticular includes the most recent calibration of WR line 
luminosities used to synthesize the WR features, up-to- 
date stellar tracks, and recent stellar atmospheres for O 
and WR stars complemented by Kurucz models for cooler 
stars (see SV98 for a full description). In all cases the 
high-mass loss stellar tracks of Meynet et al. (1994), which 
reproduce a large number of properties of individual WR 
stars and WR populations in nearby galaxies (Maeder & 
Meynet 1994) are used. Note that except for the improved 
O star atmospheres used by SV98 the recent Starburst99 
synthesis models (Leitherer et al. 1999) use the same basic 
input physics. 

The basic model parameters we consider are: 

a) Metallicity. Two sets of stellar tracks with different 
metallicities are explored here: Z = 0.02 (solar) and 
0.04. 

b) IMF slope and upper mass cut-off (M^p). We adopt 
a Salpeter IMF (slope a=2.35), and A/up=120 Mq as 
our standard model. 

c) Star formation history (SFH). Models for instanta- 
neous bursts (coeval population), extended burst du- 
rations (constant SF during period Ai; in this case 
age=0 is defined at the onset of SF, i.e. corresponds 
to that of the oldest stars present), and constant SF 
are considered. Models of combined stellar populations 
(arbitrary ages, and relative weights) with the same 
SFH have also been calculated. 

d) Fraction of ionizing Lyman continuum photons (f-y). 
fry indicates the fraction of ionizing photons absorbed 
by the gas. Our standard value is /.y = 1. Values /-y < 1 
are used to simulate various effects (e.g. dust absorp- 
tion, photon leakage outside regions, etc.) leading to a 
reduction of photons available for photoionization. 

e) Stellar extinction is kept free parameter for the 
comparison of the SED. This is justified by numer- 
ous indications showing different extinction between 
gas and stars (e.g. Fanelli et al. 1988). Note that 
such a difference affects also EW(Ha), EW(HP), and 
/(WR)//(H/3) ratios. In Figs. 4, 5, 8 we therefore show 
two sets of observations: corrected assuming C(H/3)v, = 
C(H/3)gas (filled symbols) and using C{'Rf3)^, (open 
symbols) determined from the SED model fit (see be- 
low). 

One of our main aims is to constrain the IMF (slope, 
Mup). To do so we adopt in this section a "conservative 
approach" which consists of exploring as far as possible 
a "standard" (Salpeter) IMF to see if any observational 
constraint requires a deviation from it. Regarding the WR 
features this in particular implies that in this section (i.e. 
Figs. ^ and ^ all observations will be compared to syn- 
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log EW(HjS) log EW(H/;) 

Fig. 4. Observed and predicted equivalent width (left panel) and line intensity with respect to H/3 (right panel) 
as a function of EW(Ill3). Observed points are: Mkn 309 (triangle), Mkn 589, Mkn 1199, Mkn 1259 (circles), two 
measurements for Mkn 710 (=NGC 3049) from SCK99 and GITOO (squares). Filled symbols assume a unique extinction 
correction for gas and stars; open symbols include differential stellar-gaseous extinction (see text). Data are uncorrected 
for aperture effects (cf. GITOO). Typical uncertainties are: 5-10 % for EWiRP), < 10 % for i;W^(WRbump), and - 
0.05 dex in log(/(WR)//(II/3)). Model predictions are shown for instantaneous bursts at Z=0.02 (sohd hue) and 
Z=0.04 (dashed line), and extended bursts at Z=0.02 (dotted line; burst durations At = 2,4,6,8,10,12 Myr). All 
models assume a Salpeter IMF with Mup=120 Mq. The arrow illustrates the shift which has to be applied to the 
model predictions in case of a decrease of H/3 e.g. due to absorption of photons by dust (see text). The observed WR 
bump is well reproduced by extended bursts with durations At ~ 4-10 Myr. 



can thus be taken as a lower limit for the predicted WR 
strengths which increase with increasing metallicity. The 
use of this lower limit is necessary in view of the lack of 
precise metallicity determinations for metal-rich objects. 

5.2. Results 

The comparison of the observed and predicted WR fea- 
tures is shown in Fig. |[ We show predictions for instanta- 
neous bursts at Z=0.02 (solid line), Z—O.OA (dashed), and 
burst durations of At = 2,4,6,8,10,12 Myr at Z=0.02 
(dotted). All models are calculated with the "standard" 
IMF (Salpeter, Afup=120 Mq). Figure| shows the follow- 
ing: 

— both the observed WR equivalent width and the 
/(WR) / I{il(3) ratio fall below the predictions for 
an instantaneous burst at solar metallicity rather than 
between the Z=0.02 and 0.04 curves corresponding to 
super solar metallicities. If each observed region cor- 
responded indeed to a single coeval population, this 
finding would be a strong argument in favor of a de- 



21-25 Mq) with respect to the standard IMF. How- 
ever, in the present case this can be explained quite 
simply as follows. 

— The observations are well reproduced by Z=0.02 mod- 
els assuming extended burst durations of At ^ 4-10 
Myr. The corresponding ages of our objects, as indi- 
cated by EW{ll/3), are between ~ 7 and 15 Myr as 
shown in Fig. ^. The observations of the red WR bump 
(C IV A5808; see Fig. ||) and the remaining observa- 
tional constraints are well matched by this scenario 
(see below). 

— Including the differential extinction correction for gas 
and stars as derived from the Balmer decrement and 
the SED (cf. below) leads to non-negligible changes of 
the considered quantities. Allowing for this effect does, 
however, not modify the above conclusions for the bulk 
of the objects. 

The extended burst scenario at ages ~ 7-15 Myr also 
explains the presence of red (super)giants revealed by the 
TiO bands (Sect. 4). Even at a quantitative level we note, 
somewhat surprisingly, that the broad TiO features (es- 

.i^iu, £;ocn X f^„4- i-U^ — ^"u: — j-^ . 



12 



Daniel Schaerer et al.: Massive star populations and the IMF in metal-rich starbursts 




log EW(HjS) 

Fig. 5. Same as Fig. ^ for the equivalent width of the 
C IV A5808 WR bump. 

quite well reproduced by our SEDs, where Kurucz models 
are used to describe cool stars (see Fig. 0). In any case, we 
find that the overall shape of the SED alone already indi- 
cates a mixture of populations with different ages. Inde- 
pendent of the amount of extinction, single stellar popula- 
tions are not able to produce the required relative amount 
of blue 5100 A) to red Hght. 

Last, but not least, the overall SEDs in the optical 
range of all objects are very well reproduced with the ex- 
tended burst models at Z=0.02 0. The SEDs, shown in 
Fig. 0, correspond to the model for At=8 Myr, taken at 
the age deduced from the observed H/3 equivalent width 
(see Fig. Both the predicted stellar and nebular contin- 
uum emission are included in the models. In the present 
case the latter is of importance only for emission short- 
ward of the Balmer break. However, in this domain a 
quantitative comparison is not possible since the H and 
He continuous emission coefficients depend quite strongly 
on the electron temperature which is poorly known for the 
present objects. 

As discussed above, the extinction suffered by the 
stellar continuum (C(H/3)^) is taken as a free parame- 
ter, since it may differ from the value derived from the 
gas (Balmer decrement). The resulting values of C(H/3)j, 
are systematically smaller than the values from the gas 
(Mkn 309: C(H/3)* = 0.15 / C(H/3)gas =0.68, Mkn 589: 
0.5/0.72; Mkn 710: 0.25/0.52, Mkn 1199: 0.2/0.55, Mkn 
1259: 0.2/0.60) with differences similar to those found e.g. 
by Calzetti (1997) and Mas-Hesse & Kunth (1999). For 
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Fig. 6. Temporal evolution of the predicted H/3 equivalent 
width for the solar metallicity (Z=0.02) models shown in 
Figs. ||and ||. Sohd hue: instantaneous burst, dotted lines: 
extended bursts with durations At = 2, 4, 6, 8, 10, 12 Myr. 
For Z=0.04 the same equivalent widths are obtained at 
ages on average ~ 1 Myr younger. The shaded band shows 
the range of the observed equivalent widths. 



Mkn 309, with the largest AC(H/3), we have independent 
evidence showing that C{ll(3)^, < C(H/3)gas: the observed 
ratio EW(Ha)/ EW(HP)= 8.9 clearly exceeds the max- 
imum (typical) value of ~ 7.5 (~ 5-6) predicted for in- 
stantaneous (extended) burst models^. Any other effect 
(underlying old population, multiple components etc.) will 
further reduce the predicted EW (Ha) / EW (HP) ratio; the 
only way to bring both in agreement is by adopting a 
larger extinction for the gas. Accounting for the differen- 
tial extinction and adopting the SEDs shown in Fig. 0, the 
observed EW{Ila)/ EW{}if3) are in excellent agreement 
with the model predictions for all objects. 

We conclude that all the given observational con- 
straints can be reproduced by models with a Salpeter IMF 
extending to high masses for a burst scenario with star for- 
mation extending over ~ 4-10 Myr. Let us now address 
the question if this model is unique and what other models 
can be excluded. 



® Although less clearly the remaining objects seem to indi- 
cate the same finding. EW (Ra) / EW (Rfi) is larger than the 
typical value predicted for extended burst models. Mkn 589: 
EW{'lia)/EW{}i/3)=6.30; Mkn 710: 6.68; Mknll99: 6.84; Mkn 
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5.3. Alternate models 

A variety of other models have been considered to explain 
the observations. We here consider mostly three varia- 
tions, all of them adopting a large upper mass cut-ofF: 
1) the effect of internal dust absorbing the ionizing pho- 
tons, 2) variations of the star formation history, and 3) 
variations of the IMF slope. 

The first one is illustrated in Fig. |[ where the arrow 
indicates the direction along which the predicted lines will 
be shifted if 50 % of the Lyman continuum photons are lost 
(e.g. absorbed by dust inside the H ii region). The Figure 
shows that in this case the observations can still reason- 
ably be explained with our "standard model" ; somewhat 
longer burst durations are then required. A distinction be- 
tween the two cases from the remaining constraints includ- 
ing the SED is not possible. However, values of ^ 0.5 
seem more difficult to reconcile with Fig. ^. 

As mentioned earlier instantaneous burst models are 
excluded from the observations due to the presence of both 
WR and RSG features, and the detailed SED. This also 
holds for synthesis models allowing for the formation of 
WR stars through Roche lobe overflow in massive close 
binary systems (Cervifio 1998, Cerviho et al. 2000). At the 
epoch where such instantaneous burst models predict the 
simultaneous presence of WR (formed in interacting sys- 
tems) and RSG, the predicted H/3 equivalent width is too 
low (Cerviho, private communication). Constant star for- 
mation (CSF) predicts a maximum /(WR)//(H/3) which 
is much too small (see Fig. ^), and too large H/3 equivalent 
widths. The only way reconcile this scenario is to invoke 
effect 2) with variable values of as low as ~ 0.1-0.25. 
The SED does not exclude this possibility. We have also 
performed Monte-Carlo simulations combining arbitrarily 
two "standard IMF" instantaneous bursts of different ages 
and relative contributions. Several possible combinations 
are found satisfying all the observational constraints; all 
of them obviously involve a young (^ 5 Myr) population 
superposed to a older one (^ 7-9 Myr) producing enough 
evolved red stars. Again, no strong distinction is found 
( "by eye" ) for the SEDs compared to the extended burst 
scenario discussed above. 

Finally we have considered the case of a steeper IMF 
slope of a 3.3, keeping Miip=120 M©. This IMF 
deficient in massive stars with respect to the Salpeter 
slope produces just a strong enough WR bump (up to 
£;H^(4650) ~ 9 (13) A for Z=0.02 (0.04)) at the maxi- 
mum during an instantaneous burst. Extended bursts are 
essentially excluded in this case, and very few solutions 
(combined populations) can be found. Not enough mas- 
sive stars (WR stars with Mini ~ 20-25 Mq) are expected 
for this steep IMF slope. 

From the models considered in this Section we con- 
clude that a variety of models with a "standard IMF" 
(Salpeter slope and large Mup) reproduce the observa- 



an IMF slope a = 3.3 close to the Miller & Scalo (1979) 
value; this renders such a steep IMF in metal-rich regions 
very unlikely. 

Before we proceed to a further discussion of the IMF 
we shall briefly compare our successful models with those 
available by other authors for Mkn 710 and Mkn 1259. 

5. ^. Comparison with other work 

The WR region in Mkn 710 has recently been studied in 
detail by SCK99 and Mas-Hesse & Kunth (1999). Our re- 
sults are in agreement with the study of SCK99, whose 
measurements have been included in the present paper 
(see also SCK99 for comparisons with the earlier data of 
Vacca & Conti 1992). This region being the only metal- 
rich in the sample of SCK99, no detailed fitting attempt 
was made. Regarding the IMF, SCK99 conclude that the 
Salpeter slope produces enough WR stars to explain the 
observations. From the analysis of 5 metal-poor regions 
SCK99 found that short burst durations (Ai < 2-4 Myr) 
were required; Mkn 710 was not considered. The some- 
what longer burst durations found for the present metal- 
rich sample could be due to an increasing complexity of 
the regions considered here (nuclear SB versus isolated 
H II regions) hosting progressively more separate popula- 
tions. 

For the comparison with Mas-Hesse & Kunth (1999) 
it is important to note that their optical-UV observations 
include a much larger region (10 by 20") than the present 
observations typically extracted over 2-3" wide regions. 
This explains most likely the larger contribution of an 
older underlying population contributing to the optical 
flux found by these authors. The results on the burst du- 
ration are in good agreement. Incidentally Mas-Hesse & 
Kunth derive the same stellar extinction as found from 
our work. 

For Mkn 1259 a subsolar metallicity was derived by 
Ohyama et al. (1997) based on [O ill] A5007A. The value 
used in the present work, based on a spectrum of superior 
S/N and resolution by GITOO, should be more accurate. 
Using their low metallicity value, Ohyama et al. (1997) 
concluded that an instantaneous burst with an age of ~ 5 
Myr reproduces well the observed properties of the nuclear 
region of this object. Using our metallicity, the strengths 
of the WR features imply an extended burst, compatible 
with the observed RSG features and the SED. For most 
aspects however, the difference with Ohyama et al. (1997) 
is probably negligible. 

6. Discussion 

Having excluded steep IMF's with a high Mup, we will now 
examine what limits can be set up on Mup for a Salpeter 
slope and discuss several implications from our results. We 
then briefly discuss the plausibility of the star formation 
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Fig. 8. Predicted dependence of the equivalent width of 
the WR bump on variations of the upper mass cut-off from 
Mup= 120 Mq (thin hnes) to Mup = 30 M© (thick Unes). 
Instantaneous burst models with a Salpeter IMF for so- 
lar metallicity (dotted, solid) and Z=0.04 (short-dashed, 
long-dashed) are shown. The shaded area indicates the 
maximum domain covered by twice solar metallicity mod- 
els with the reduced mass cut-off (Mup = 30 Mq). Upper 
panel: Temporal evolution of EW{WIV). Lower panel: Pre- 
dicted EW (WR) as a function of EWIJilS). The observa- 
tions are shown using the same symbols as in Fig. IJ. The 
strong dependence on Mup illustrates the power of using 
the WR features to constrain Mud- Discussion in text. 



6.1. Constraints on the IMF from the present data 

What value of Mup is our data exactly sensitive to ? From 
the ages and burst durations derived above we see that 
the youngest stars present have ages of the order of ~ 3-5 
Myr. This corresponds to stars with initial masses Mini ^ 
40-120 Mq (see Meynet et al. 1994). Although in the best 
case our data could be sensitive to changes of Mup over 
this range, we have to conservatively admit that we might 
not be able to distinguish between different values for the 
upper mass cut-off ^ 40 Mq . To improve the present situa- 
tion one obviously has to be able to study the intrinsically 
youngest regions and/or include additional observational 
constraints. 

With the objects studied here it is not straightforward 
to answer the reverse question, namely "What lower limit 
can we set on Mup ?". In the ideal case (instantaneous 
burst with a known metallicity) the use of WR lines is a 
very powerful discriminant, as illustrated in Fig. 8. Here, 

;„ ™„j-„n;„;j-,. j- / v n r\A\ „i u 



be adopted to obtain a conservative lower limit for Mup. 
However, as it is clear from Sect. 5, the use of the WR 
features is not straightforward for answering the above 
question given the extended duration of SF. On the other 
hand, as apparent from Fig. H, the observed H/? equivalent 
widths indicate clearly an upper mass cut-off of Mup > 30 

Mq. 

The /(He i A5876)//(H/3) ratio has previously been 
used by Bresolin et al. (1999) as an indicator of Mup. We 
note that for all of the data except Mkn 309 our limit 
on Mup seems also consistent with the observed /(He i 
A5876)//(H/?) ratio as estimated from simple Case B re- 
combination theory. More firm conclusions await com- 
bined stellar population and photoionization modeling (cf. 
Sect. 7). 

6.2. IMF constraints from Ha/H(3 equivalent widths of 
young SF regions 

Assuming a Salpeter IMF with Mup as large as ~ 100-120 
Mq at high metallicity, to be compatible with the obser- 
vations discussed above, has obviously implications which 
need to be confronted with other observational data. A 
straightforward implication is that at age=0 H ii regions 
are expected to show very large H/3 equivalent width ('^ 
450- 500 A for Z=0.02-0.04, and larger values for Z < 
Zq and Mup=120 Mq). It is well known that regions 
with such large EW{Il(3) are apparently not found or 
are extremely rare, independently of metallicity to first 
order. What is the situation for high metallicities ? Rel- 
atively few EW{Il(3) measurements have been published 
for metal-rich H ii regions. Among the data of Shields et 
al. (1991), Oey & Kennicutt (1993), Kennicutt & Garnett 
(1996), van Zee et al. (1998), and Bresolin et al. (1999, 
see Bresolin 1997) few regions with 12 -I- log(0/H) ^ 8.9 
have EW{}i/3) ^ 250 A. Furthermore these measurements 
are most likely affected by very large uncertainties due to 
placement of the weak continuum as inspection of some 
data of van Zee et al. shows. From this data we conclude 
that ma,x{EW(Hf3)) ~ 200-300 A for metal-rich H ii re- 
gions. For a Salpeter slope this corresponds to Mup ~ 35 
- 50 Mq , compatible with our conservative lower limit for 
Mup but clearly below "normal" values of ^ 100-120 Mq. 

Although a detailed discussion and interpretation of 
the observed Ha or H/3 equivalent widths distributions is 
beyond the scope of the present work, we would like to 
briefly comment on four processes which could lead to re- 
duced values of the observed EW(H(3) maximum: 1) Dust 
inside the H II regions, 2) underlying populations diluting 
the continuum, 3) early phases hidden in compact radio 
sources (cf. Kobulnicky & Johnson 1999), 4) absence of a 
proper ZAMS for high mass stars in the accretion scenario 
of Bernasconi & Maeder (1996). From the available data 
we are not able to discuss possible selection effects. 

The first two are the most commonly considered pos- 
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objects, remains to be explored. For metal-poor objects 
(e.g. H II galaxies) with a low dust content the effect 
should on average be small. Although for metal-poor ob- 
jects the young bursts might dominate the visible light 
and underlying absorption (2) should thus not be impor- 
tant for the max(£'H^(H/3))-discrepancy (cf. Mas-Hesse & 
Kunth 1999), the recent study of Raimann et al. (2000) 
suggests the existence of non-negligible underlying popu- 
lations. If the bulk of star formation (SF) occurred in a 
similar fashion as in He 2-10, where Kobulnicky & John- 
son (1999) find a significant number of compact radio 
sources, this could constitute a possibility to hide the earli- 
est phases of SF and thereby reducing the observed max- 
imum EW{llf3). We note, however, that the estimated 
lifetime of these regions (~ 10 % of the O star lifetime; 
Kobulnicky & Johnson) is clearly too short to resolve the 
max(i?VF(H/3))-discrepancy which concerns the first ~ 2 
Myr. Effect 4 has been suggested by Mas-Hesse & Kunth 
(1999) as a possibility to explain this discrepancy. How- 
ever, test calculations (Schaerer & Bernasconi, 1995, un- 
published) show that the stellar tracks following the ac- 
cretion scenario of Bernasconi & Maeder (1996) lead to an 
insufficient decrease of max (SVF(H/3))Q In addition, the 
validity of their scenario remains to be verified by detailed 
confrontations with observations of individual stars. 

Regarding only metal-rich region we conclude from the 
available EW{llf3) data that the observed maximum H/3 
equivalent width is compatible with a Salpeter IMF with 
an upper mass cut-off Mup ~ 35 - 50 Mq, in agreement 
with the results from our WR galaxies. Larger values of 
Mup cannot be excluded from our present (small) sample 
of objects, in contrast to the indications from EW{llf3) 
measurements in larger samples. In view of the discrep- 
ancy between the observed and predicted max{EW (Jif3)) 
which seems to hold for all metallicities and for which cur- 
rently no clear explanation exists, interpretations based on 
the "exclusive" use of EW{}i(3) should, however, be taken 
with caution. 

6.3. Are the derived burst durations plausible ? 

An interesting result derived from our models is that all 
the objects analysed here require extended burst durations 
over ~ 4-10 Myr or a superposition of two or more bursts 
with similar ages differences and a young component (^ 
5 Myr). ^ For the most distant object from our sample 
(Mkn 309) the size of the observed region corresponds to 
~ 2 kpc. For a typical galaxian velocity dispersion of 200 
km s^^ the crossing time is 8.4 Myr. Shorter times are 
obtained for the remaining objects (cf. Heckman et al. 

^ max(£;W(H/3)) ~ 400 A for Z=0.02 and the standard IMF. 
We note that the observed H/3 equivalent width and the 
/(WR)//(H/3) ratio of IRAS 01003-2238 (cf. refs. in Table |) 
can be explained by Z—0.04 models with a longer burst (At ^ 
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1997 for Mkn 477) including also IRAS 01003-2238, which 
has a very compact morphology (Surace et al. 1999). This 
simple estimate shows that the above burst durations pose 
no fundamental physical problem for any of the objects 
discussed here. 

Is there any difference of increasing burst duration 
with distance as could be expected on average if the ob- 
served regions are composed of an increasing number of 
"small" super star cluster like objects ? In this context it is 
interesting to note that in an earlier study of WR galaxies 
with distances of 4-14 Mpc burst durations not exceeding 
~ 2-4 Myr were found using the same evolutionary syn- 
thesis models (SCK99). This difference is likely due to two 
effects: first the distance effect, and second the different 
nature of the objects included in the two samples (dwarf 
galaxies with one or few SF regions versus predominantly 
nuclear starbursts). 

7. Summary and conclusions 

We have obtained new spectroscopic observation of Mkn 
309, a starburst galaxy with one of the largest WR popu- 
lation known (cf. Osterbrock & Cohen 1982, Armus et al. 
1988). Our high S/N observations allow in particular a de- 
tailed determination of the WR star content of this object, 
composed of WNL, WCL and possibly also some WCE 
stars (Sect. 4). Adopting the extinction from the Balmer 
decrement (but cf. below) this corresponds to massive star 
populations of - 25000 WNL, ~ 10000-15000 WC, and 
^ 65000 equivalent 07V stars in an area of ^ 1.7 x 3.5 
kpc. Several methods to estimate the metallicity concur 
to indicate 12 -I- log(0/H) ~ 9.3-9.4, which makes Mkn 
309 the most metal-rich WR galaxy known today. 

The Mkn 309 observations have been combined with 
spectra of four other metal-rich WR galaxies from the 
sample of Guseva et al. (2000) to constrain the basic prop- 
erties of the massive star populations (IMF slope, Mup) 
and the star formation history (age, burst duration) of 
these objects by detailed quantitative comparisons with 
appropriate evolutionary synthesis models (Sect. 5). We 
have used the following main observational constraints: H 
recombination lines, intensities and equivalent widths of 
the main WR features, TiO bands originating from a pop- 
ulation of late type stars (which are most likely red super- 
giants; cf. Sect. 4) and the detailed SED from - 3500-7000 
A. Considering a large parameter space for the synthesis 
models (metallicity, IMF slope and upper mass cut-off, 
star formation history, possible effects of dust or photon 
leakage on the ionized gas, variable stellar extinction) we 
obtain the following main results: 

1.) The observations are explained by continuous star 
formation extending over ~ 4-10 Myr ("extended bursts") 
seen at ages of 7-15 Myr or a superposition of several 
bursts with similar age differences (4-10 Myr) including a 
young (< 5 Myr burst). This naturally explains both the 
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and the presence of red supergiants. Such extended burst 
durations, compared to other WR galaxies indicating even 
shorter SF timescales (e.g. SCK99), are plausible in view 
of the physical sizes of the observed regions and the nature 
and morphology of our objects (nuclear starbursts), and 
pose no fundamental physical problem (Sect. 6). 

2. ) The SEDs in the optical range are very well repro- 
duced for all objects, provided the stellar light suffers from 
a smaller extinction than that of the gas (derived from the 
Balmer decrement). Differences of AE{B — V) = 0.15-0.36 
are found for our objects with a mean Ay = 1.3±0.1. Our 
finding of smaller stellar extinction based solely on optical 
data is in agreement with earlier studies (e.g. Fanelli et al. 
1988, Calzetti 1997, Mas-Hesse & Kunth 1999) requiring, 
however, the combination of UV and optical data. 

3. ) The observational constraints are compatible with 
a Salpeter IMF extending to masses AI^p 40 Mq . Anal- 
ysis of younger regions are required to provide real upper 
limit on the value of the upper mass cut-off. We have also 
explored what deviations from this "universal" IMF are al- 
lowed. Adopting a conservative approach we derive a lower 
limit of Mup~ 30 Mq for the Salpeter IMF. Larger val- 
ues (Mup^ 30-40 Mq) are obtained using probably more 
realistic assumptions on the metallicity and SF history. 
An upper mass cut-off of at least ~ 35 - 50 Mq seems 
also indicated from H/3 equivalent width measurements 
of metal-rich H ii regions in spiral galaxies (see Sect. 6). 
Quite independently of the value of Mup, steeper IMFs 
with slopes such as the Miller & Scalo (1979) IMF, are 
very unlikely. 

Obviously the solutions fitting the observable con- 
straints are not unique. The main conclusions just sum- 
marised are, however, quite robust to changes of the con- 
sidered model parameters (see Sect. 5^). Nevertheless it 
should be recalled that the results presented here rely 
quite extensively on our present knowledge of massive 
stars, especially on the evolution of WR stars. The ad- 
vantage of studying objects of solar or somewhat higher 
metallicity is that this corresponds to the domain with the 
largest number of individual WR stars in the Local Group 
which has previously been compared quite successfully to 
the evolutionary models used here (see Maeder & Meynet 
1994). In this sense the application of our synthesis mod- 
els to starburst galaxies should at least allow meaningful 
differential comparisons of properties such as the IMF. 

Through the analysis of the massive star content in lu- 
minous metal-rich starbursts by means of their direct stel- 
lar features, our aim is to constrain the upper end of the 
IMF (slope and upper mass cut-off) in metal-rich environ- 
ments in general and in particular in massive starbursts 
such as ULIRGs. Whereas previous studies based on indi- 
rect indications from the ionized gas seem to indicate an 
IMF deficient in massive stars (e.g. Goldader et al. 1997, 
Luhman et al. 1998, Bresolin et al. 1999, but see the recent 



do not support such important differences compared to 
the "standard" Salpeter IMF. 

To obtain more reliable constraints on the IMF in 
metal-rich environments new high quality spectral obser- 
vations will be necessary. While analysis of few individual 
objects are extremely useful, larger samples are e.g. re- 
quired to constrain the exact value of M^p by sampling all 
relevant ages. To include the combined information from 
gas and stars, coupled evolutionary synthesis and pho- 
toionization models including all the stellar and nebular 
observables should be used consistently. Furthermore the 
statistical distributions of the properties of e.g. metal-rich 
disk H II regions also contain encoded information about 
the IMF, SF history etc. The theoretical tools for such 
analysis are available today. With the new powerful tele- 
scopes available now, great progress should be made on 
studies of stellar populations and the IMF in metal-rich 
environments. 
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Fig. 7. Comparison of observed spectral energy distributions with predicted ones for extended burst models (A<:=8 
Myr) with a standard IMF. All SEDs are normalised at 4800 A. The predicted SEDs, which do not include emission 
lines, are reddened by the variable amount C(H/3)^, indicated in each panel. 



